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The spin polarization P of the transport current through tlie interface between superconducting 
Al and ferromagnetic Fe is determined by means of Andreev reflection at nanostructured point 
contacts. We observe a systematic decrease of P with decreasing contact resistance. Our data 
provide evidence for the reduction of P by spin-orbit scattering and thus establish a link between 
density-of-states and transport spin polarizations. 
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The spin-dependent electronic transport in magnetic 
heterostructures is an emergent field because of poten- 
tial applications in spintronics [l], Q • The knowledge of 
the current spin polarization P of materials is a key issue 
for the functionality This is particularly important 
for high-P materials such as half-metallic alloys, e.g., 
Heusler compounds [1, 0] . The spin polarization of fer- 
romagncts can be measured by various techniques, in- 
cluding spin-polarized photoemission , spin-dcpcndcnt 
electron tunneling 01, and point-contact Andreev reflec- 
tion (AR) between a ferromagnet F and a superconduc- 
tor S 0] . An issue of considerable importance is how the 
current spin polarization obtained by AR l9| is related to 
the ferromagnct's bulk spin polarization [lO|. 

AR has been widely used to extract the current s pin 
polarization in a variety of materials [III. [isl . [T^ . Iisl . 
|l6[. In particular, Upadhyay et al. [I4I employed the 
type of contacts also used by us. AR at an interface be- 
tween a normal metal N and S is the coherent process by 
which an electron from N incident on S is retro-reflected 
as a hole, thus creating a Cooper pair in S. The theoreti- 
cal analysis of these experiments is usually carried out in 
the spirit of the Blonder-Tinkham-Klapwijk (BTK) the- 
ory [131 ■ The ordinary reflection at the N/S interface is 
parametrized by a phenomenological parameter Z. The 
sensitivity of the Andreev process to the spin of the car- 
riers originates from the spin content of a Cooper pair 
(5 = in conventional superconductors) and the conser- 
vation of the spin direction at the interface. Employing 
a ferromagnet as the nonsupcrconducting electrode leads 
to a reduction of the AR probability for spin-singlet su- 
perconductors because incoming electron and outgoing 
hole must have opposite spins Q. Two major approaches 
have been suggested to extract P from AR datas: (i) the 
AR current is decomposed into a fully spin-polarized and 
fully unpolarized current fldl . [isl . [l9j , (ii) spin-dependent 



transmission coefficients rt and T| are introduced, hence 
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P =1 T| — T| I /{ti +Ti) [20] ■ For a comprehensive recent 
review of different models to extract the AR spin polar- 
ization, see Rcf.im In this letter, we report on AR spec- 
tra obtained on a series of nanostructured Al/Fe contacts 
that show a distinct dependence on contact size. A possi- 
ble origin of this systematic size dependence is spin-orbit 
scattering, shedding light on the origin of the differences 
in the spin polarization determined by different experi- 
mental methods. We chose Fe in this work because it 
does not show surface degradation of P when prepared 
under UHV, in contrast to many Heusler compounds that 
show a reduced P at the surface. 

We employ nanostructured Al/Fe point contacts (PCs) 
that are formed by two metal films connected by a tiny 
hole in an insulating membrane in between [2^ . These 
contacts are very stable and reproducible spectra are ob- 
tained even after several thermal cycles. Furthermore, in- 
situ preparation of the metal films in a single run avoids 
oxidation or degradation of the interface between the two 
electrodes. Holes of radius a between 7.5 and 25 nm were 
fabricated by e-beam structuring of a 50- nm thick Si3N4. 
a was determined from the contact resistance (see below) 
after evaporation in UHV of a 200-nm thick Al layer on 
the side of the membrane with a bowl-like funnel towards 
the hole, and of a 12-nm thick Fe layer topped by a Cu 
layer of thickness dcu = 188 nm on the opposite (flat) side 
(see Refs.llO, [22| for details). By variation of the process 
parameters for the e-beam lithographical prestructuring 
and subsequent reactive ion-etching, contacts with resis- 
tance between 1 and « 30 17 are obtained. Low- 
ohmic PCs {R^ < 2^1) implying a large contact diame- 
ter usually show clear evidence of non-ballistic transport 
and therefore were not used for the present investigation. 
The differential conductance G ~ dl/dV vs. applied bias 
voltage V was measured in a dilution refrigerator between 
T « 20 mK and 2K. 

Fig.[T] exemplarily shows the differential conductance 
G{V) of a PC with i?Ar(4.2K) = ISAQ at various tem- 
peratures T below Tc of the Al film. G is normalized 
to the normal-state conductance Gm ~ ^n^- The data 



* Electronic address: 
t Electronic address: 



gemot . goUQphys . uni-karlsruhe . de] 



h. vlQphys . uni-karlsruhe. de| 



2 




J I I I I I I 1 L 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

V(mV) 

FIG. 1: (Color online) Differential conductance G = dl/AV 
normalized to the normal-state value Gn = G{T > Tc) — 
54.3 mS for a nanostructured Al/Fe point contact (sample 
#5) for different T. For clarity, the curves are shifted upwards 
successively by 0.2 units with decreasing T. The conductance 
curves calculated with the T|- T| model match the experimen- 
tal curves exactly except for the lowest T in the region of the 
G maximum (solid line). 



are analyzed in terms of two spin-dependent transmis- 
sion coefficients T| and for the majority and minority 
charge carriers in F, respectively. We assume that the 
PC can be described with a single pair of transmission 
coefficients [l^, H^l which contain the microscopic prop- 
erties relevant for the transport, i. e., the spin-split band 
structure of F, the electronic structure of S, and the in- 
terface properties. For t < 1, an electron experiences 
ordinary reflection with a finite probability 1 — r. This 
leads to a minimum of G at zero bias and sharp maxima 
at voltages corresponding to the superconducting energy 
gap A(r). Using T|, r^, and A as free parameters for 
each T, the model yields excellent least-squares fits to 
the Andreev spectra over the whole T range below Tc 
and the entire voltage range, as shown in Fig.[TJ Only 
at the lowest nominal T = 17mK, where the fit yields 
T| = 0.348, rx = 0.997, and A = 166 fxeV, can the fits 
be distinguished from the data, the calculated maxima 
appear somewhat sharper than experimentally observed. 
This might be caused by a levelling-off of the electron T 
due to heating by electromagnetic stray fields or a small 
pair-breaking effect by Fe. The excellent agreement be- 
tween experiment and theory is also reflected in the T de- 
pendence of the zero-bias resistance Gq^ = {<il / dV)Y^Q 
shown in Fig.jSK- P is found to be 0.482 for this partic- 
ular sample. The characteristic parameters for a total of 
six contacts are listed in TableHl 

Shown in Fig. [5] are G{V) curves obtained at T w 0.1 K 
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FIG. 2: (Color online) Normalized differential conductance 
G/Gn of nanostructured Al/Fe point contacts at T ~ O.IK. 
For clarity, the curves are shifted upwards successively by 0.15 
units. The conductance curves (solid lines) are calculated 
with the r|- Tj model (see text). 



TABLE I: Superconductive energy gap Ao = A(r ^ 0), tran- 
sition temperature T^^'^ , transmission coefficients Tf , r;, and 
current spin polarization P as determined by a fit of the An- 
dreev spectra at T < 200 mK for point contacts with contact 
resistance Rm ~ 1/Gn at T = 4.2 K. a is the contact radius 
calculated from the Wexler resistance (see text). 
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Rn 


a 


Ao 




r| 


n 


P 


No. 


m 


(nm) 


(meV) 


(mK) 








1 


2.68 


24.7 


0.174 


971 


0.371 


0.983 


0.452 


2 


6.98 


14.3 


0.175 


903 


0.362 


0.984 


0.460 


3 


7.29 


13.9 


0.157 


1009 


0.349 


0.993 


0.480 


4 


9.59 


12.0 


0.190 


970 


0.361 


0.984 


0.462 


5 


18.4 


8.4 


0.166 


1174 


0.348 


0.997 


0.482 


6 


24.2 


7.3 


0.174 


1115 


0.343 


0.994 


0.487 



for all samples. One clearly sees a systematic reduction 
of the amplitude of the AR-related features when going 
from sample 1 (small Rn) to sample 6 (large Rn), in- 
dicating a systematic increase of P. Again, the quality of 
the fits is excellent over the whole voltage range up to 
V = 4Ao/e. 

The high stability of our PCs, together with the high 
quality of the fits, allows a detailed investigation of the 
T dependence of A(T) and P{T) as displayed in Fig.jSh 
and [3h. T is normalized to the experimental of 
each sample in order to account for the minor variation 
of Tc from sample to sample. T^'^p has been determined 
from the sudden drop of the zero-bias resistance Gq^{T) 
(see Fig. 13^). A(T) follows the standard BCS dependence 
(solid line in Fig.[3}3). At first sight, this might be surpris- 
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FIG. 3: (Color online) (a): T dependence of the zero-bias re- 
sistance Gq ^ together with a theoretical curve (see text), (b): 
T dependence of the superconducting gap A for all samples 
investigated in this study, obtained from a least-squares fit of 
the data and normalized to A(r < 200 mK < Tc). (c): T 
dependence of the current spin polarization P ~\ — \ 
/{ti + Ti). T is normalized to the experimental T^^'^ of each 
sample in order to account for the small variation of Tc from 
sample to sample. 



ing because the superconductor is in direct contact with 
a ferromagnet. Apparently, pair breaking plays a minor 
role as also inferred from the small rounding of the AR 
spectra mentioned above. Here, the small dimensions of 
the S/F contact are instrumental, because they limit the 
AR laterally to a tiny region and thus do not strongly af- 
fect S. Indeed, Tc appears to be somewhat lower for low 
i?N, i-e., larger contacts (see TableU]). The current spin 
polarization P for each sample is independent of T within 
the error of the least-squares fitting procedure (Fig. [31;). 
Only close to Tc where the spectra become very shallow 
does the scatter in P increase. The finding that P is 
independent of T strongly supports our assignment of P 
as an intrinsic parameter of the particular S/F PC under 
study. 

Inspection of Table 1 reveals that indeed P depends on 



Rn in a systematic fashion, i.e., larger contact resistances 
go along with larger P. We interpret this finding as fol- 
lows. In view of the fact that all samples were prepared 
in the same fashion, we are confident that it is chiefly 
the contact radius a that determines Rn. In particu- 
lar, we recall that the interface between Al and Fe was 
never exposed to air. The radius a of a circular homo- 
contact can be calculated from i?jv for different regimes 
of current flow through the contact, depending on the 
ratio of a to the electron mean free path i of the ma- 
terial, with the limiting cases of the ballistic or Sharvin 
regime if a <^ ^ and the thermal or Maxwell regime if 
a ^ £. For the intermediate regime an interpolation 
for Rn = "f{K)p/{2a) + 4pif/(37ra) was suggested by 
Wexler [23] ■ Here ^{K) is a slowly varying function of 
K = l/a with 7(0) = 1 and -i{K oo) = 0.694 and 
the Maxwell {K — > 0) and Sharvin limits {K — * oo) are 
given by R'^-^ = p/2a and R§^°° = Api/Sna^ fH, 
respectively. Using arithmetically averaged parameters 
pi = 2.66 • 10-15 flrn^ and p = 5.84 • lO'^ flm [2^, and 
solving the Wexler formula iteratively for self-consistent 
values of K, a is found to vary between 7.3 and 24.7 nm, 
i.e., the change in Rn corresponds to a change in a by a 
factor of 3 (see Table 1). 

P{a) decreases systematically as displayed in Fig.Sl 
The error bars reflect the statistical error of typically 6-8 
measurements, all taken below 200 mK w 0.2 for each 
sample. We note that the overall trend of P{a) docs not 
depend on the details of how a was calculated. A pos- 
sible explanation of the P{a) dependence is as follows. 
For large a the incident electrons probe a larger volume 
before being retro-reflected as a hole than for a small a. 
Any scattering process affecting P will therefore have a 
stronger effect for larger contacts [l^- We suggest that 
spin-orbit scattering — with a constant scattering length 
— is operative as the main source of the reduction of P, 
and model the dependence P{a) by a simple exponential 
decay, P{a) = Pq exp (— a/^so)- We obtain a spin-orbit 
scattering length £so = 255±91nm from a fit to the data 
in Fig.m This value is very reasonable albeit somewhat 
lower than the spin diffusion length Xf = 500 - 1000 nm 
found for Al [28j , which ~ apart from the error introduced 
by the scatter of the data - may be attributed to a smaller 
elastic mean free path of our Al films or minor interdif- 
fusion at the Al/Fe interface [2§|. We note that using 
the model of Mazin et al. [l3 with the fit parameters A, 
P, and Z, the same systematic trend of P(a) with even 
smaller Iso is obtained with no systematic correlation be- 
tween P and Z . With this interpretation, our data consti- 
tute the first instance of a clear relation between the spin 
polarization of the Andreev current and the bulk spin po- 
larization. Pq = 0.496±0.01 almost reaches the value ob- 
tained by spin-polarized photoemission spectroscopy on 
polycrystalline Fe films [30| and on the (111) surface of 
single-crystalline Fe [Slj yielding P = 0.54 and P = 0.6, 
respectively. The domain structure of the Fe film is ex- 
pected to play no appreciable role for our PCs because 
both the domain size and the thickness of domain walls 



4 



0.50 




a (nm) 

FIG. 4: Dependence of the current spin polarization on the 
contact radius. The soUd hue shows P = Po exp {—a/£so) 
with Po = 0.496 ±0.01 and £,o = 255 ±91 nm. Over the short 
distance a <^ £so the exponential decay of P{a) looks almost 
like a straight line. 



are much larger than the contact size |32l | . 

For Pb/Co and Pb/Ni nanostructured PCs the current 
polarization was reported to be independent of the con- 
tact size with AG/Gn = ±0.2 at any (fixed) voltage p^ . 
We note that our systematic variation of AG/Gn = ±0.2 



at any voltage is also of that order. However, given the 
fact that Pb exhibits larger spin-orbit scattering, a larger 
variation would be expected. Possibly, the different be- 
havior originates from a different microscopic structure of 
the interface and/or disorder in the PC region, perhaps 
from intcrfacial spin-dependent scattering whose impor- 
tance was demonstrated by investigations of the influ- 
ence of the Co film thickness on the spin-filtering effect 
in nanostructered Pb/Co contacts [33|. Indeed P ob- 
tained for the Pb/Co and Pb/Ni contacts is substantially 
lower than P of our Al/Co [201 and Al/Fe contacts. This 
underscores the need to prepare samples under identical 
conditions in order to determine the dependence of the 
current spin polarization on the contact radius as done 
in the present work. 

In conclusion, we have found a systematic dependence 
of the spin polarization P of the Andreev current in S /F 
contacts on the contact size. We suggest that this size de- 
pendence arises from spin-orbit scattering which reduces 
P and is most effective in large contacts. This size ef- 
fect must be taken into account when determining P via 
Andreev reflection. 

We thank J. C. Cucvas. W. Bclzig, and K. Samwcr for 
helpful discussions. 
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